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ABSTRACT 


GOES  West  Infrared  (IR)  and  Visual  (VIS)  satellite  observations  for 
the  North  Pacific  Ocean  area,  3S-55N,  120W-1800,  at  2345  GMT  for  eight 
June/July  1978  dates  are  processed  to  yield  20  statistical  parameters 
which  are  analyzed  for  their  use  in  discerning  the  existence  of  marine 
fog.  The  exploratory  sample  of  522  satellite  observations , at  3 nmi  x 
3 nmi  resolution  (at  subsatellite  point)  .are  related  to  the  associated 
0000  GMT  synoptic  ship  reports  which  serve  the  role  of  ground  truth.  The 
best  discrimination  between  fog/nc-fog  observations , using  the  statisti- 
cal parameters,  occurs  for  stratification  of  the  data  by  the  meridional 
component  of  the  wind  (south  vs  north),  IR  temperature  (>  269  K vs 
< 269  K)  and  the  standard  deviation  of  IR  temperature  (~  2 K vs  > 2 K) . 
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I.  INTRODUCTION  AND  BACKGROUND 


Ship  accidents  due  to  reduced  visibility  because  of  fog  exact  a heavy 
toll  in  materiel  and  human  life  (Quayle,  1976).  For  the  Navy,  marine  fog 
curtails  carrier  fli^rt  operations  as  well  as  surface  craft  operations. 

A study  over  a five-year  period  concerning  dollar  cost  and,  more  impor- 
tantly losses  of  human  life,  emphasizes  the  impact  of  marine  fog  upon 
naval  operations  (jiheeler,  13.74). 

Accurate  diagnosis  and  prognosis  of  .marine  fog  regions  over  the  open 
■ocean  would  reduce  the  type  of  casualties  referred  to  above  and  enhance 
the  capability  of  any  naval  operation  in  which  the  limiting  parameter  is 
horizontal  visibility.  A specific  application  of  fog  prognosis  is  in 
Optimum-Track  Ship  Routing  (3TSR) , a program  in  existence  at  Fleet  Numeri- 
cal weather  Central,  Monterey,  California  CfTJWC). 

The  only  operational  objective  hemispheric  fog  forecasts  generated 
on  a daily  basis  are  those  produced  by  the  FNWC  for  all  areas  of  the 
Northern  Hemisphere.  Their  product,  called  FTER,  is  a computerized  pro- 
bability of  advective  fog,  and  is  promulgated  twice  daily  at  0000  and 
1200  GMT  for  forecast  intervals  up  to  72  hours.  These  fog  forecasts  are 
statistically  derived  from  numerical  model  output  parameters  (U.S.  Naval 
Weather  Service  Cormvand , 1975). 

Observations  by  meteorological  satellites  may  be  ’useful  in  Improving 
the  diagnosis  of  marine  fog  areas  and  concomitant  increases  in  forecast 
accuracy.  Research  conducted  at  the  Naval  Postgraduate  School,  Monterey, 
California  (NPS),  in  1974  and  1975  , utilized  visual  (VIS)  and  infrared  (IR) 
weather  satellite  imagery  from  polar-orbiting  satellites  to  qualitatively 
diagnose  the  presence  of  marine  fog  over  the  North  Pacific  Ocean.  Initial 
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results  indicated  that  further  investigation  into  the  ’use  of  digital 
satellite  data  for  determining  marine  fog  was  necessary  (Vial lace  and 
Renard,  1975).  Using  MOM- 2 digital  satellite  data,  continuing  investi- 
gations in  1975  established  optimal  ranges  of  digital  count  values  (DC/'s) 
or  gray  shades,  for  fog  occurrence  for  both  VIS  and  IR  data  (Hale  and 
Renard,  1975).  The  level  of  ski'-  achieved  in  applying  these  findings 
indicated  that  further  research  was  necessary.  Research  conducted  in 
1975  and  1977  was  directed  to  the  feasibility  of  using  VIS  and  IR. 
digital  data  from  both  the  military  polar-orbiting  meteorological  satel- 
lite, Defense  Meteorological  Satellite  Program  (DMSP)  Satellite,  and  the 
Synchronous  Meteorological  Satellite  C3MS-2 ) , as  recorded  over  the 
eastern  Perth  Pacific  Ocean  area  tlhli  and  Renard,  1977).  Optimal  VIS 
and  IR  ranges  for  the  occurrence  cf  marine  fog  were  determined  for  SMS-2 , 
while  the  IMS?  data  set  was  too  limited  for  meaningful  analysis.  The 
results  obtained  in  this  study  were  promising  enough  to  warrant  a more 
extensive  test  -using  Geostationary  Operational  Environmental  Satellite 
(GOES)  data.  Such  a study  was  initiated  in  fall  1973  utilizing  GOES  West 
June/ July  1973  digital  data  from  the  eastern  Morth  Pacific  Ocean. 


— 


OBJECTIVES  ART  APPROACH 


The  immediate  objective  of  the  subject  research  is  to  identify  physi- 
cal statistical  par am; tens  derivable  from  VIS  and  IR  weather  satellite 
observations  which  cay  be  used  to  discriminate  ocean  areas  of  marine  fog 
from  areas  where  no  fog  exists.  The  ultimate  objective  is  to  combine 
such  parameters  with  numerically-derived  model  output  parameters  in  a 
regression  scheme  to  specify  fog  areas  at  initial  time  and  yield  a per- 
sistence factor  useful  in  prognostic  time.  The  numerical  statistical 
asoect  af  the  research  has  been  on-gcing  since  1377  (Van  3 man  and  Pen.ard , 
1977;  Quinn,  1973). 

The  approach  utilized  to  pursue  the  immediate  objective  was  to  study 
various  statistical  distributicns  cf  parameters  dependent  cn  temperature 
ar.d  reflected  radiance  as  observed  by  radiometric  sensing  in  the  I?,  and 


713  channels,  respectively. 


h.e  hcrth  Jaci:ic  Ice  an  ICET-Vest  satel- 


lite. Related  statistical  approaches  were  utilized  by  the  In  ter- American 
Tropical  Tuna  icmmissicr.  (Stevenson,  Kirkhara  and  Madsen . 12”7)  to  screen 
clouds  for  sea-surface  temperature  diagnosis  and  by  the  university  of 
Maryland  (Tooth,  13” 3 ) for  objective  clcud  type  classification.  These 
research  endeavors  provided  important  background  for  this  study. 

Tata  from  polar-orbiting  satellites,  such  as  SUSP,  were  net  ccnsiderec 
here,  in  part  due  to  the  experience  of  the  crevicus  researchers  (Ihli  and 
Renard,  1977).  The  GOES  has  the  advantage  of  taking  observations  at 
short  intervals  and  hence  can  be  readily  matched  in  time  to  ground  truth 
synoptic  obser/atiens . 


limn 
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A.  .AREA  OF  STUDY 

The  area  of  study  was  limited  to  that  portion  of  the  eastern  Morth 
Pacific  Ocean  within  the  following  boundaries:  35-55N;  120W-170E 
(Figure  1).  The  area  selected  was  a compromise  due  to  several  factors: 
the  coverage  provided  by  the  GQES-West  satellite;  data  storage  constraints 
imposed  by  the  Space  Science  and  Engineering  Center,  University  of 
Wisconsin,  Madison,  Wisconsin  CSSEC),  where  the  data  were  received;  and 
the  climato logical  area  of  marine  fog  occurrence  (Figure  1). 

3.  PERIOD  OF  STUDY  AND  GROUND-TRUTH  DATA 

July,  climato logically  speaking,  records  the  most  frequent  occurrence 
of  marine  fog  over  the  Morth  Pacific  Ocean  (U.S.  Navy,  Naval  Oceanography 
and  Meteorology,  1978).  Thus  it  was  arranged  to  have  GOES  West  observa- 
tions recorded  at  the  SSEC  for  the  period  starting  with  0245  GMT  27  June 
1973  and  ending  with  2345  GMT  4 August  1973.  Specifically,  VIS  and  IR 
digital  data  were  retrieved  daily  for  0245,  1745,  2045  and  2345  GMT  with 
certain  exceptions  (.Table  I).  The  exceptions  were  due  in  part  to  the 
navigational  parameters  being  unavailable  for  all  times  during  this 
period,  and  machine  and  operator  error.  The  amount  of  data  actually 
used  was  a subset  of  3 days  between  27  June  and  13  July  1973.  This 
period  was  selected  since  the  operational  satellite  was  transferred  to 
a stand-by  mode  while  replacing  it  with  a new  vehicle  in  the  14-16  July 
time  period.  Thus,  it  appeared  that  a sufficient  exploratory  sample  was 
available  for  study.  The  possibility  of  greater  variance  in  the  satel- 
lite parameters  due  to  combining  data  from  two  satellites  was  thus  avoided. 
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The  0000  GMT  primary  synoptic  observation  time  only  was  selected  for 


study  since  all  longitudes  under  study  were  in  daylight  at  that  time 
(1000  LST  at  120W  to  1100  LST  at  165E)  and  VIS  as  well  as  IR  data  (at 
2140  GMT)  could  be  considered.  'Hie  surface  ship  weather  observations 
which  provided  "ground  truth"  were  obtained  from  FNWC.  These  observa- 
tions were  classified  into  one  of  41  fog  likelihood  categories  ( FLO ) 
(following  Willms,  1975)  which  were  further  combined  into  seven  TLC 
Groups  (FLCG)  and  filially  three  divisions  (.FLCGD)  (Tables  IT,  III). 

This  was  accomplished  for  the  1001  ship  reports.  After  considering  all 
of  the  fog  cases  (FLGG  1-0)  (Figure  7)  „uui  selecting  only  those  corre- 
sponding non-fog  cases  (TLCG?)  in  the  s.une  general  vicinity,  the  number 
of  ship  reports  was  reduced  to  522,  or  approximately  5 2 A of  the  original 


data. 

C.  GOES  WEST  DATA 

GOES  West  digital  data  were  furnished  on  magnetic  tape  by  the  TEC 
in  their  save  tape  format , a means  ot  recording  manv  data  frames  .into  a 
single  .magnet ic  tape.  The  dynamic  range  of  both  the  VIS  and  TR  sensors 
is  0 to  255  (NASA,  1979).  lew  (high)  DCV  values  for  VIS  data  correspond 
to  surfaces  having  lew  thigh)  albedo,  while  low  (high)  DCV  values  for  TR 
data  indicate  relatively  high  tlcw)  temperatures. 

Tie  Visible  and  Infrared  Spin  Scan  Radiometer  (VISSR)  aboard  GOES 
West  scans  in  the  west/east  direction  on  each  spin  of  the  spacecraft. 

In  the  visible  region  (9. 55-0. "5  micrometers)  VISSR  scans  eight  parallel 
lines  on  each  spin,  providing  about  1 km  (0.5  nmi)  resolution,  while  in 
the  infrared  (10.5-12.6  micrometers)  VISSR  scans  one  line  on  each  spin, 
providing  roughly  4 km  (4  nmi)  resolution  at  nadir-  (MOM,  NASA,  19"6). 
These  line  to  line  resolutions  deteriorate  u;  the  field  of  view  (FOV) 
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moves  further  away  from  the  subpoint  (figure  3).  Along  each  scan  line 
the  data  resolution  is  nominally  1 km  (.0.5  nmi)  for  VIS  data  and  4 km 
(2  nmi)  for  IR  data  at  the  subsatellite  point  (equator  at  135W) , with 
less  deterioration  than  the  line  to  line  resolution  at  points  removed 
from  nadir. 

Since  ship  positions  are  recorded  to  tenths  of  degrees  and  since 
areas  of  marine  fog  were  deemed  to  be  quite  extensive,  it  was  decided 
that  a 3 nmi  x 3 nmi  pixel  size  would  be  sufficient  for  the  processing 
of  the  digital  satellite  data  thereby  reducing  the  quantity  of  data  from 
that  available  and  minimizing  the  ensuing  cost  of  processing.  The  3 nmi 
VIS  data  are  obtained  by  averaging  successively  every  six  samples  (each 
sample  at  half  nmi  intervals)  along  a line,  placing  the  result  at  the 
fourth  sample  position.  In  the  case  of  IR  where  every  other  line  is 
redundant  information,  every  third  element  and  line  are  moved  by  one  mile 
where  elements  and  lines  are  at  2 nmi  intervals  (Figure  4). 
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IV.  DATA  ANALYSIS  AND  RESULTS 


A.  APPROACH 

The  approach  adopted  was  to  compute  20  statistical  parameters  for 
both  the  VIS  and  IR  observations  (Appendix  A,  after  3ooth,  1973)  at  and 
in  the  vicinity  of  the  locations  of  the  522  0000  GMT  ground-truth  synop- 
tic ship  reports  and  to  compare  the  average  values  for  each  of  seven 
FLCG's.  These  7 PLUG’S  were  combined  into  3 divisions  (Tables  II,  III) 
in  order  to  have  a minimal  number  of  related  groups,  thereby  allowing  a 
sufficient  number  of  ship  observations  to  be  included  in  each  division 
so  as  to  give  statistical  meaning  to  the  parameter  values. 


3.  GOES -WEST  DATA  ANALYSIS 


Satellite  orbital  elements  necessary  to  earth  locate  the  satellite 
observations  were  provided  by  the  SSEC.  These  were  received  in  final 
usable  form  at  the  NPS  in  early  January  1979.  Prior  to  this  time,  during 
a five-day  period  at  the  SSEC  in  August  lS^S,  five  days  of  the  grciand 
truth  synoptic  ship  report  data  were  located  in  satellite  coordinates 
(line,  pixel)  by  the  author.  Tie  process  of  retrieving  line  and  element 
positions  for  each  ship  observation  position  on  the  Me IDAS  Interactive 
Uompu  or  at  the  SSEC  proved  to  be  inefficient  for  the  author's  purposes 
.and  only  provided  the  location  of  the  ship  report  in  each  data  frame  time 
(i.e.,  time  of  receipt  of  satellite  ooservations) . This  process  dis- 
allowed further  efficient  use  of  the  VIS  and  IR  data  sets  for  each  frame 
time  and  thus  did  not  contribute  to  the  computation  of  the  20  statistical 
parameters.  Moreover,  it  was  discovered  in  late  UY78  that  an  incomplete 
set  of  navigational  parameters  was  available  from  the  SSEC.  Fortunately, 
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the  satellite  division  at  the  Naval  Environmental  Prediction  Research  Fa- 
cility, Monterey,  California  OIEPRF)  was  provided  with  the  3SEC  software 
for  earth  locating  the  satellite  data,  thus  it  became  more  expedient  to 
perform  the  navigational  processing  at  the  NPS  for  all  of  the  satellite 
data.  Ship  position  (latitude  and  longitude),  total  cloud  cover,  low  cloud 
cover  and  type , height  of  low  cloud  base,  middle  and  high  cloud  types, 
surface  wind  direction  and  speed,  FLC,  and  FLCG,  along  with  the  satellite 
frame  time,  were  coded  as  inputs  to  the  analysis  program. 

Both  the  7 IS  and  IR  GOES-West  digital  satellite  data  were  checked  cn 
each  of  the  eight  days  for  the  distribution  of  DCV 's  in  the  range  0-255. 

Due  to  data  processing  at  )IES3,  "holes"  or  cerocounts  appeared  methodically 
in  the  distribution  of  IR  DCV's  (Figure  5).  However,  this  problem  was 
rectified  by  converting  DCV-values  into  equivalent  black-body  temperatures 
according  to 


XV 

330  - ~~  , XV  < 1?5 


TOO  = 


413  - DCV  , DC/  > 176 


This  conversion  was  dene  before  statistically  manipulating  the  data  to 
derive  the  20  statistical  parameters.  This  problem  was  net  exhibited  by 
the  VIS  data. 

Further , in  the  IR  data  a temperature  of  33C  K occasionally  occurred. 
Such  a daren  was  considered  an  indicator  of  missing  or  erroneous  data, 
thus  invalidating  the  'usefulness  of  the  particular  .array  of  satellite 
data  surrounding  a ground- truth  point.  .Also,  since  the  temperature  range 
300-330  K has  no  meteorological  significance  in  the  area  of  study,  IR  DC/ 
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values  between  0 and  6Q  are  suspect.  A similar  statement  is  true  for 
the  temperature  range  118-163  K (JR  D CV  values  between  220  and  255). 
Therefore  any  non-zero  DCV  in  either  of  the  above  ranges  may  be  attribut- 
able to  signal  processing  noise  of  the  GOES-West  satellite  digital  data. 

Next  the  ship  position  was  located  inside  the  data  frame  in  both  the 
VIS  and  IR  data.  Initially,  an  area  approximately  84  x 90  nmi  at  45M, 

160W  (IS  lines  x 32  pixels)  surrounding  the  ship  location  was  extracted 
from  the  data  frame,  and  the  statistical  parameters  computed  for  both 
IR  and  VIS  data.  This  area  size  appeared  to  be  too  large,  especially 
for  a FOV  near  55N,  160W  Ci.e. , 119  x 105  nmi).  The  area  size  was  reduced 
to  3 lines  by  16  pixels,  which  is  approximately  42  x 45  nmi  at  45N,  150W, 
and  60  x 52  nmi  at  55N.  160W.  This  decision  was  not  entirely  arbitrary 
since  a compromise  between  an  area  of  sufficient  size  to  justify  statis- 
tical averaging  and  yet  small  enough  to  be  representative  of  events  at 
the  ground- truth  point  was  necessary.  Considered  to  be  a partial  justi- 
fication for  the  size  selected  is  the  fact  that  the  averages  of  the  array 
mean  temperatures  and  of  the  array  mean  albedos  were  not  very  different 
from  the  averages  of  the  temperatures  and  of  the  albedos  located  at  the 
ship  position. 

Having  computed  the  2Q  parameters  for  each  observation  point,  scatter 
diagrams  were  used  to  compare  the  VIS  and  IR  parameters  in  order  to  deter- 
mine seme  discrimination  between  fog  (FDCG1)  and  no-fog  (FLCG7).  There 
appeared  to  be  little  success  here.  (Per  example  compare  Figure  6(a)  and 
6(b).)  These  early  efforts  were  largely  motivated  by  the  need  to 
separate  cases  of  lew  cloud  only  (with  or  without  fog)  from  cases  "con- 
taminated" by  middle  and/or  high  cloud.  With  this  in  mind  a stratifica- 
tion of  oases  according  to  average  IR  temperature  cut-offs  of  260  K and 


17 


269  K was  investigated.  The  260  K cut-off  was  abandoned  after  shewing 
little  evidence  of  being  able  to  achieve  discrimination  of  fog/no  fog 
cases.  The  269  K cut-off  appeared  to  be  a more  likely  discriminator, 
at  least  in  principle,  since  700  kPa  climatological  temperatures  indi- 
cate this  to  be  a reasonable  low  value  across  the  study  area  of  the 
North  Pacific  Ocean  (.Figure  7),  and  thus  a suitable  minimum  for  tops  of 
low  clouds.  Moreover,  the  distributions  of  fog/ no  fog  cases  from  the 
limited  sample  are  nearly  equal. 

Next  the  269  K mean  IR  temperature  was  combined  with  the  standard 
deviation  of  temperature,  S ^ 2 K,  as  a trial  discriminator  of  fog/no 
fog  cases  (Figure  3).  Thus,  statistics  for  the  following  array  mean 
temperature  and  standard  deviation  CT,  3m)  categories  were  tabulated  for 
the  ground- truth  synoptic  ship  report  locations:  all  T;  T > 269  K; 

T > 269  K with  S-  < 2 K;  T > 269  K with  3m  > 2 K-,  T < 269  K;  T < 269  .< 
with  3m  2 K;  and  T < 269  K with  3m  > 2 K (Table  IV) . At  this  point 
there  appears  to  be  little  discrimination  between  fog/ no  fog  divisions 
I,  II,  III  for  each  temperature/ standard  deviation  stratif ication . Of 
course  there  are  differences  in  the  values  of  the  statistical  parameters 
from  one  temperature/ standard  deviation  category  to  another  but  this  is 
to  a large  measure  merely  evidence  of  expected  variation  due  tc  differ- 
ences between  middle/high-cloud  contaminated  cases  (T  < 269  K)  and 
mostly  low-cloud  arrays  (T  _>  259  K) , and  in  part  due  to  variable  cloud 
type  in  each  stratification  (as  stratified  vs  cumuliform  clouds). 

The  next  step  was  to  stratify  the  sample  by  cloud  cover  (M),  a vari- 
able contained  in  the  synoptic  report  at  each  ground-truth  location.  It 
was  reasoned  that  some  veiling  of  discrimination  between  fog/no  fog  was 
due  to  variable  sky  cover  and  hence  variable  T,  3,  gradient,  etc.  Simply 
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the  reports  were  separated  by  N >_  .8  and  M < .3  with  statistics  shown 
for  the  former  (Table  V).  Unfortunately,  the  sample  is  considerably 
reduced,  but  this  is  mostly  due  to  elimination  of  no -fog  cases  and  the 
attendant  partially  cloudy  skies;  the  effect  is  considered  positive.  In 
this  test  there  appears  to  be  increased  fog  discrimination  among  the 
parameters,  more  specifically  between  FLCG1  and  FLCG7.  However,  the 
sample  for  T < 263  K with  S^,  < 2 K is  too  small  to  be  considered  signi- 
ficant for  this  analysis. 

Ihe  test  was  repeated  'using  as  a discriminator,  where  <_  M.  As 
before,  the  reports  were  separated  by  _>  .3  and  < . 3 with  statistics 
shown  for  the  former  (Table  VI). 

Previous  work  in  the  Department  of  Meteorology,  'IPS  (Retard  and 
Servaas,  1377;  Van  Drman  and  Renard,  1377)  indicated  the  value  of  wind 
direction  as  a discriminator  of  fog/ncrfog.  Thus,  the  satellite  data 
were  further  processed  to  relate  the  IR  and  VIS  samples  to  ground-truth 
observations  with  .3  and  south  (100°  to  2T°)  or  north  (230°  to 

290°)  winds  (calm  cases  were  emitted).  Ihe  best  results  to  date  were 
obtained  in  this  segment  of  the  study  (Tables  VII,  VIII).  In  order  to 
focus  on  the  increasing  parameter  discrimination  of  fog/nc  fog  cases 
with  each  enhanced  stratification  step  (Tables  V-VIII),  the  parameter 
differences  Division  I (PLCG1)  - Division  III  (FLCG7)  were  computed  and 
are  presented  in  Table  IX.  It  is  obvious  from  the  table  that  the  best 
division  of  fog  (FLCGl)/no  fog  CFLC37)  occurs  with  T _>  269  K and 
3,  < 2 K for  south  and  north  wind  segments.  For  example , without  the 

L 

aid  of  any  of  the  20  statistical  parameters,  for  T 769  K,  3,  < 2 K and 
south  winds,  no  fog  (FLCG7)  occurs  only  23%  of  the  time  while  north  winds 
are  associated  with  no  fog  62%  of  the  time.  The  results  are  similarly 


definite  for  T > 269  K,  S > 2 K and  south  winds,  but  the  sample  is  too 

1 

small  to  be  significant.  In  any  case,  the  statistical  parameter  differ- 
ences Div  I versus  Div  III  also  maximize  for  these  groups,  especially 
in  the  VIS  observation,  thus  allowing  for  further  modulation  of  these 
empirical  probabilities  based  on  a processing  of  the  differences. 

The  final  experiment  was  conducted  to  increase  the  fog  (Div  I) 
sample  by  allowing  FLCG2  to  be  combined  with  FLCG1,  since  FLCG2  observa- 
tions are,  in  general,  fog  at  observation  time,  although  light.  Table  X 
shews  the  differences  (Div  I^CFLCGl , 2 ) - Div  III  (FLCG7))  for  the  same 
categories  as  Table  LX  south  and  north  wind  discrimination  samples  only. 
Tonparing  the  values  of  Tables  IX  and  X,  the  advantage  of  combining 
FLCG1  -and  2 appears  questionable  and  this  facet  of  the  research  was  not 
pursued  further  at  this  time. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  principal  objective  of  this  study  was  to  identify  physical  sta- 
tistical parameters  derivable  from  VIS  and  IR  weather  satellite-observa- 
tions  which  may  be  ’used  to  discriminate  ocean  areas  of  marine  fog  from 
areas  where  no  fog  exists.  Through  the  interpretation  and  processing 
of  digital  GOES-West  satellite  data  using  conventional  surface  ship  re- 
ports as  ground-truth  verification,  a surface  report  possessing  a mean 
temperature  (T)  greater  than  or  equal  to  269  K,  a standard  deviation  of 
temperature  (3~)  less  than  or  equal  to  2 K,  and  a south  wind  component 
has  a relatively  high  frequency  of  fog  for  M.  i .3.  .Also  for  north  winds 
having  the  same  temperature  characteristics , no  fog  has  a relatively  high 
frequency.  Further,  the  differences  in  mean  statistical  parameter  values 
between  fog/no-fog  cases  tend  to  maximice  especially  for  VIS  observa- 
tions, given  scuth-ccmponent  winds,  T 2.  “59  K,  3^  £ 2 X,  M.  >_  .3.  To  a 
lesser  extent  this  statement  is  also  tree  for  north-component  winds , all 
other  conditions  the  same.  Beyond  this,  lack  of  a sufficiently  large 
sample  precludes  definite  conclusions. 

Though  the  sice  of  the  sample  is  small,  important  hints  as  to 
the  usefulness  of  GOES  West  in  indicating  marine  fog  areas  are  evident. 
Therefore,  the  following  recommendations  are  offered  for  future  studies: 

(1)  Conduct  a similar  study  utilicing  at  least  the  remaining  16 
days  of  2345  GMT  GOES-West  satellite  digital  data. 

(2)  Compute  kurtosis  values  for  the  VIS  DCV  and  IR  temperature 
distributions  as  a twenty-second  parameter. 

(3)  Further  examine  statistical  relationships  between  conventional 
synoptic  ship  report  data  and  both  VIS  DCV  and  IR  temperature 
distributions . 


Mu 
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(4)  Utilize  the  linear  regression  approach  to  test  the  significance 
of  those  satellite  parameters  (.predictors)  apparently  contri- 


buting to  specifying  the  probability  of  fog/ no  fog. 

(5)  Vary  the  stratification  IS  temperature  and  standard  deviation 

temperature  values  to  maximize  the  differences  in  mean  parameter 
values  between  fog/no-fog  cases . 


A 


APPENDIX  A:  PARAMETER  COMFJTATIOMS 


Each  set  of  3 x 16  IR  and  VIS  satellite  observations  centered  on  the 
ground-truth  GGCQ  GMT  synoptic  ship  report  is  treated  as  a one-dimen- 
sional array  of  123  values  in  computing  the  following  parameters  (refer 
to  Appendix  3 for  parameter  listing) . 

The  mean  is  computed  by: 


n 


m = 


X. 


Cn=123)  Parameter  1 


where  X-  is  the  observed  value. 

The  second  and  third  moments  about  the  mean  are  computed  by  the 
following  formula: 

n 

u.  = y (X.  -m)^  for  k=2S3  . 
i=l 


Values  of  the  standard  deviation  (.3),  coefficient  of  variation  (CV),  and 
skewness  GK)  are  computed  by  the  following  formulae: 


CV  = S/m 


Parameter  3 
Parameter  4 


and  SK  = Ug/Ci^  Parameter  5 

Values  at  1,  16,  50,  34  and  39%  cumulative  frequencies  CCF)  are 
determined  CParameters  6-10)  with  the  range  defined  as: 


range  CR)  - 39h 


Parameter  11 


Paraneter  12  is  the  median  minus  the  mean. 

The  next  six  parameters  (Parameters  13-13)  are  computed  by  consider- 
ing the  following  element  configuration: 

A T 

Ax 

D E — 3 •,  Ax  = 1 unit, 


where  the  gradient  C2)  and  Laplacian  Cl)  are  computed  by 

G = [J-  (3-D)2  + i CC-A)2]'^ 

and  L=A+3+CtD-4E. 

The  minimum  and  caximum  are  algebraically  selected  for  both  G and  1, 
and  the  arithmetic  mean  values  for  both  are  computed  in  the  usual  sense. 

The  power  spectrum  is  computed  by  transforming  each  set  of  3 :<  16 
observations  f(I,J)  into  the  spatial  frequency  domain  F(u,v)  by  using  a 
two-dimensional  fast  Fourier  transform  CF7T) . The  transformed  observa- 

A 

ticn  is  then  normalized  by  multiplying  by  its  complex  conjugate  F (y,v): 

A 

C(y,v)  = F (y,v)  • FCy,v)  . 

The  power  scectrum  parameters  CP)  (Parameters  19-21)  are  computed 
by  combining  row  and  column  variances  centered  about  the  two-dimensional 
plane  at  wavenumber  i: 


Legend: 
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APPENDIX  B:  LIST  OF  GOES-WEST  WEATHER  SATELLITE 
PARAMETERS  AND  THEIR  SYMBOLIZATIONS 
WITH  UNITS 

indicates  average  of  a number  of  single  array  values 


(single  array  mean  values) 


Parameters  Symbols  Units 

1.  Temperature  (.VIS  DCV)  T(A)  K(DCV) 

2.  Station  Temperature  (VIS  DCV)  T'sTA^STA''  K(DCV) 

2.  Standard  deviation  of  S„,(S  ) K(DCV) 

temperature  CVIS  DCV)  - A 


4. 

Coefficient  of  variation  of 
temperature  CVIS  DCV) 

CVTC CVA> 

None 

5 . 

Temperature 

CVIS  DC/)  skewness 

sEjCsic^) 

ft 

6. 

1%  cumulative  frequency  of 
temperature  CVIS  DCV) 

crf1%(CFA1%) 

K(DCV) 

7. 

16% 

t» 

ft 

3. 

50%  ’’ 

It 

Crr50%(CFA50%> 

rr 

3. 

34%  " 

Tt 

crra4%<CFA34%> 

r» 

10. 

99% 

11 

CFr99%(CFA99%) 

rr 

11. 

Temperature 

CVIS  DCV)  range 

RpCRA) 

IT 

12. 

Temperature 

median-mean 

CVIS  DCV) 

Mt(Ma) 

It 

13. 

Temperature 

gradient 

CVIS  DCV) 

5rCGA) 

K/2Ax 

(dcv  n Ax) 

14. 

Minimum  temperature  CVIS  DCV) 
gradient 

n (C.  ) 

,jmint^°mina; 

ft 

15. 

Maximum  temperature  CVIS  DCV) 
gradient 

^MAXr(GMAXA) 

*T 

16. 

Laplacian  CVIS  DCV)  gradient 

K(DCV) 

17. 

Algebraic  minimum  temperature 
CVIS  DCV)  Laplacian 

LMINt(  LMINa  5 

T A 

ft 

26 


Algebraic  maximum  temperature 
(VIS  DCV)  Laplacian 

K(DCV) 

Normalized  pcwer  spectrum 
variance  of  temperature 
(VIS  DCV):  wavenumber  1 

•T  UA  J 

^ of  total 
variance 

."formalized  power  spectrum 
variance  of  temperature 
(VIS  DCV):  wavenumber  2 

~2^p— 2) 

T UA 

r» 

Normalized  power  spectrum 
variance  of  temperature 

(VIS  DC V ) : wavenumber  3 

P?C*?> 

»» 

ABLE  I 


Available  GOES-West  satellite  navigational  para- 
meters as  a function  of  time  (GMT)  for  study  period. 


x indicates  gamma  values  available  for  picture  time 
indicates  gamma  values  missing  for  picture  time 
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TABLE  II.  Marine-fog  likelihood  categories  (FLC)  as  dependent  on 
the  visibility-weather  group  elements  of  the  primary 
synoptic  report. 


Fog 

Coded 

Coded 

Likelihood 

Coded 

Present 

Past 

Categories 

Visibility 

Weather 

Weather 

1 

ft 

40-49 

4 

2 

A 

40-49 

5 

3 

ft 

10 

4 

4 

ft 

11,12 

4 

5 

ft 

10 

5 

6 

ft 

11,12 

5 

n 

/ 

90-93 

AA 

4 

3 

90-93 

U 

4 

9 

94 

** 

4 

10 

94 

H 

4 

11 

90-93 

ft  ft 

5 

12 

90-93 

pr 

5 

13 

94 

ft* 

5 

14 

94 

H 

5 

15 

23 

4 

16 

* 

23 

5,6 

17 

.4 

40-49 

6 

13 

A 

10 

6 

19 

* 

11,12 

5 

20 

90-93 

AA 

5 

21 

90-93 

H 

5 

22 

94 

AA 

5 

23 

94 

a 

5 

24 

ft 

40-49 

2 

25 

ft 

10 

2 

26 

* 

11,12 

r) 

L. 

27 

90-93 

AA 

0 

tm 

23 

94 

2 

29 

A 

23 

2 

30 

95-96 

AA 

4 

31 

* 

40-49 

0,1, 3, 7-9 

32 

ft 

10 

0,1, 3,7-9 

33 

A 

11,12 

0,1,3, 7-9 

34 

90-93 

AA 

0,1, 3, 7-9 

35 

94 

AA 

0,1, 3, 7-9 

36 

ft 

23 

0,1,3,' 7-9 

37 

95-96 

A A 

0-3, 5-9 

38 

90-94(95-961 

H 

0-3, "-9(0-3, 5-9) 

39 

97-99 

**,H 

0-3, 5-9 

40 

97-99 

*4 

41 

95-96 

H 

Li 

4 

* Denotes  any  visibility  allowed,  i.e.,  90-99 

**  Denotes  any  present  weather  code  except  10,11,12,28,40-49  or  heavy 
present  weather  codes 

H Denotes  heavy  present  weather  codes  only:  30-39,59,62-65,67,59,72- 
75,31,82,84,86,38,90-99 


"n^BLE  III  Foe  likelihood  categories  (FLO,  groups  (FLCG)  and 
divisions  (FLCGD)  as  used  in  study. 


; — (FLCGD) 

Description 

FLC 

1 (I) 

Heavy  fog  likely  at  observation 
time  at  station 

1,2,4,6,7,11,17, 

19,24,26,27,31,33 

2 (II) 

Light  fog  likely  at  observation 
tine  at  station 

3,5,9,13,18,25,32 

3 (II) 

Fog  likely  nearby,  but  not  at 
observation  time /station 

15,16,29,36 

4 (II) 

F’robably  near  fog  in  space  and/ 
or  time 

3,10,30,40,41 

5 (II) 

Chance  of  fog  at  station  .and/or 
nearby  station 

20  , 22  , 2S, 34, 35 

6 CII) 

Visibility  qualifies  for  very 
light/ligfrt  fog,  but  not 
reported 

37 

7 (III) 

Fog  very  unlikely 

12,14,21,23,38,39 

30 


cases 


TABLE  17(a).  Statistical  parameters  for  IR  digital  data,  no  cloud  cover 
restrictions  (see  Appendices  A and  B for  parameter  informa- 
tion) . 
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TABLE  IV (b).  Continued  from  IV (a). 
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TABLE  r/(c).  Statistical  parameters  for  VIS  digital  data,  no  cloud 
cover  restrictions. 
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TABLE  IV(d) . Continued  from  IV (c). 


All  T >_  269K  T > 269K  _ T < 269K  T < 269K 
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TA3LE  V(a).  Statistical  parameters  for  IP,  digital  data,  for  M >_  .3. 
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CABLE  VII (c).  Statistical  parameters  for  VIS  digital  data  for  M and 

M.  >_  .8,  south-component  wind  directions  (100°  to  270°). 
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TABLE  VII(d).  Continued  from  VII(c). 
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TABLE  VIII(a).  Statistical  parameters  for  IR  digital  data  for  M and 
•>  .3,  north-component  winds  (230°  to  090°). 
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TABLE  V IIKb).  Continued  from  VIIKa) 
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Figure  2.  Locations  of  248  ('round-truth  observations  for  FI0TL  through 
(0000  (m  2B , 20  June,  1,  2,  3,  0,  0,  14  .July  107B). 
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gure  4.  Scheme  for  selection  of  full  resolution  IR  and  VIS  GOES-West 
satellite  data  to  fulfill  3 nmi  x 3 nmi  data  criterion.  For 
VIS,  (Rvalues  used;  for  IR  O and  Q values  used. 
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Frequency  distribution  of  TR  DCV  distribution  for  data  frame 
time  234f>  GMT  27  June  1J78.  DCV  of  zero  (25S)  is  represented 
by  upper  left  (lower  right)  comer  square . 
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Figure  5(a).  Comparison  of  average  IF  temperature  vs  average  VIS  CCJ  for 
fog" observations  (F1C1),  no  temperature  restriction. 
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Figure  5(b).  Comparison  of  average  IR  temperature  vs  average  VIS  CCV  for 
no-fog  CFLC7)  observations,  no  temperature  restriction. 
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Figure  3(a).  Comparison  of  average  IR  temperature  vs  standard  deviation  o 
temperatures  for  fog  observations  (FLC1),  restricted  to  T > 


4-i  C i 


JR  TEMP.  STD.  DEV. 

U.  00  2.00  tj.00  6.00  8.00  10.00  12.00 


i to 


LIST  OF  REFERENCES 


Booth,  A.  L. , 1973:  Objective  Cloud  Type  Classification  Using  Visual 
and  Infrared  Satellite  Data.  Technical  Mote  BN  763,  Institute  for 
Fluid  Dynamics  and  Applied  Mathematics,  University  of  Maryland, 

62  pp. 

Hale,  R.  E.  and  Renard,  R.  J. , 1975:  Use  of  the  NOAA-2  Digitized 
Satellite  Data  for  Diagnosing  Marine  Fog  in  the  North  Pacific 
Ocean  .Area.  NPS-51H175091,  Naval  Postgraduate  School,  Monterey, 
California,  92  pp. 

Dili,  C.  3.,  Jr.  and  Renard,  R.  J. , 1977:  The  Use  of  DMS?  and  SMS-2 
Digital  Satellite  Data  for  Identifying  Marine  Fog  in  the  Eastern 
North  Pacific  Ocean  .Area.  NPS-63Di77C31,  Naval  Postgraduate  School, 
Monterey,  California,  99  pp. 

National  Oceanic  and  Atmospheric  Administration  and  National  Aeronautics 
and  Space  Administration,  1976:  The  30ES/SMS  User's  Guide.  Editors: 
R.  ?.  Ccrbeli,  0.  J.  Callahan  and  W.  J.  Kotsch , NCAA/ NASA, 
Washington,  D.  0.,  113  pp. 

Ouayle,  R.  G. , 1976:  leather  and  Maritime  Casualty  Statistics. 

Mar,  lea.  Log,  Vol.  20,  No.  2,  pp.  74-80. 

fulnn , ?.  F. , 1973:  Further  Development  of  a Statistical  Diagnostic 
Model  of  Marine  Fog  Using  FNWC  Model  Output  Parameters.  M.  S. 

Thesis  Naval  Postgraduate  School,  Monterey,  California,  73  op. 

Renard,  R.  J.  and  Servaas,  T. , 19’77:  Climatological  Fog  Frequencies 
as  a Function  of  Wind  Direction  and  Speed,  Unpublished  research, 
Department  of  Meteorology,  Naval  Postgraduate  School,  Monterey, 
California. 

Stevenson,  M.  R. , Kirkham,  R.  G.  and  Madsen , 3.  J.,  19  : Development 

and  Testing  of  a Cloud  Screening  Technique  for  Use  with  Satellite- 
home  Scanning  Radiometers.  MERER  Cent  .Nc.NCC'2  3-"S-C-313  3 , Inter- 
American  Tuna  Com.,  Scripts  Inst,  of  Ccg.  , La  folia,  CA.,'95  pp~~ 

U.  S.  Navy,  Naval  Oceanography  and  Meteorology,  1978:  Study  of  World- 
wide Occurrence  of  Fog,  Thunderstorms , Supercooled  Low  Clouds  and 
Freezing  Temperatures  (MAVAIR  5G-1C-60CHI) . NWSED,  Federal  3uil±Lng, 
Asheville,  N.C. 

U.  3.  Navy,  Naval  Weather  Service  Command,  19 ”0:  Selected  Level 

Heights,  Temperatures  and  Dew  Points  for  the  Northern  Hemisphere 
(MAVAIR  50-1C-52).  Superintendent  of  Documents,  U.S.  Government 
Printing  Office,  Washington,  D.C. 


65 


U.  3.  Navy,  Naval  Weather  Service  Command,  1975:  Numerical  Environmental 
Products  Manual , MA.VAIR  50-1G-522. 

Van  Orman,  3.  L.  and  Renard,  R.  J. , 1977:  Statistical  Diagnostic  Model- 
ling of  Marine  Fog  Using  Model  Output  Parameters.  M.  S.  Thesis, 
Department  of  Meteorology,  Naval  Postgraduate  School,  Monterey, 
California,  94  pp. 

Wallace,  R.  T.  and  Renard  R.  J. , 1975:  The  Use  of  Meteorological  Satel- 
lites for  Discerning  Marine  Fog.  NPS-51Wa75031,  Naval  Postgraduate 
School,  Monterey,  California,  59  pp. 

Wheeler , 3.  E.  and  Leipper,  D.  F.,  1974;  Marine  Fog  Impact  on  Naval 

'Operations.  MPS  Report  53Wh74091,  Department  of  Oceanography,  Naval 
Post  graduate  School,  Monterey,  California,  118  pp. 

Willms,  G.  P. , 1975:  A Climatology  of  Marine-Fog  Frequencies  for  the 
North  Pacific  Ocean  Summer  Fog  Season.  M.  S.  Thesis,  Department  of 
Meteorology,  Naval  Postgraduate  School,  Monterey,  California,  59  pp. 


INITIAL  DISTRIBUTION 


Defense  Documentation  Center 
Cameron  Station 
Alexandria,  Virginia  22  314 

Library,  Cede  0142 

Naval  Postgraduate  School 

Monterey,  'California  93940 

Department  Chairman , Code  53 

C\3. ry£  ^ f ro  Lc 

Naval  Postgraduate  School 
Monterey  t.-a ",  \ fornix  339^0 

Dr.  Robert  I.  Renard,  Code  53Rd 
Department  of  Meteorology 
Naval  Postgraduate  School 
Monterev , California  33343 

LI  Otto  P.  Me Nab 
1325  Mehle  Avenue 
.Arab!,  Louisiana  ~G032 

Commanding  Officer 

Fleet  Numerical  Weather  Central 

Monterev , 339^0 

Commanding  Officer 

Naval  Environmental  Prediction  Research  Facilitv 
Monterey , California  33940 

Department  of  Meteorology,  Code  53 
Naval  Postgraduate  School 
Monterey , California  3 3943 

Department  of  Oceanography,  Code  53 
Naval  Postgraduate  School 
Monterey , California  3 3940 

Director,  Naval  Oceanography  and  Meteorology 
National  Space  Technology  Laboratories 
Bay  St.  Louis,  Mississippi  39520 

Mr.  Ihcmas  0.  Haig 

Space  Science  and  Engineering  Center 

University  of  Wisconsin 

1225  West  Dayton  Street 

Madison,  Wisconsin  5 3 ”06 


No.  Copies 


o 


LT  Carl  3.  Dili,  USN 

U.3.  Fleet  Weather  Car crax , 3uam 

3ox  12,  COMMAVMARIANAS 

FPO  San  Francisco  96630 

Dr.  Glenn  H.  Jung,  Code  63 
Department  of  Oceanography 
Naval  Postgraduate  School 
Monterey,  California  93940 

Dr.  Dale  F.  Leipper,  Code  68 
Department  of  Oceanography 
Naval  Postgraduate  School 
Monterey,  California  93940 

Mr.  Paul  R.  Lowe 

Naval  Environmental  Prediction  Research  Facility 
Monterey,  California  93940 

Dr.  Alan  Weinstein 

Naval  Environmental  Prediction  Research  Facility 
Monterey , California  9.3940 


